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Recent experiments on the S=l/2 antiferromagnetic chain compound, Cu benzoate, discovered 
an unexpected gap scaling as approximately the 2/3 power of an applied magnetic field. A theory of 
this gap, based on an effective staggered field, orthogonal to the applied uniform field, resulting from 
a staggered gyromagnetic tensor and a Dzyaloshinskii-Moriya interaction, leading to a sine-Gordon 
quantum field theory, has been developed. Here we discuss many aspects of this subject in consider- 
able detail, including a review of the S=f/2 chain in a uniform field, a spin-wave theory analysis of 
the uniform plus staggered field problem, exact amplitudes for the scaling of gap, staggered suscep- 
tibility and staggered magnetization with field or temperature, intensities of soliton and breather 
peaks in the structure function and field and temperature dependence of the total susceptibility. 

I. INTRODUCTION 

The effect of a magnetic field on an S=l/2 antiferromagnetic chain has been extensively investigated theoretically 
over many years. The Hamiltonian is written: 

H = J2[ J $i ■ ~ 9VbHS*}. (1.1) 

3 

An important conclusion was that the groundstate remains gapless right up to the saturation field. The low-energy 
excitations can be described by bosonization which predicts gapless excitations at wave-vectors and 7r and also at 
the incommensurate wave-vectors 

jfei = ±2irm{H) and k 2 = tt ± 2irm{H), (1.2) 

where m(H) is the magnetization per site, m =< Sf >. The first, detailed experimental study of such systems at 
large fields with g^sH of 0{J) were only performed very recently,!!! on Cu benzoate. This material has a relatively 
small exchange energy, J 1.57 meV so that gusH/J w .52 for a field of 7 T. While these experiments verified, in 
detail, the expected field-dependent shift of the wave-vector at which a gap minimum occurs, they also discovered an 
unexpected result. A non-zero gap Appeared which seemed to scale as approximately H 2 / 3 , with strong dependence on 
the field orientation. Dender et al.u suggested that this gap might arise from the staggered g (gyromagnetic)-tensor, 
associated with the low symmetry of the crysta l str ucture and the presence of 2 crystallographically inequivalent Cu 
sites on each chain. Thus the last term in Eq. (|l . l| ) must be replaced by: 

Hh = ~Hb £ H s [g u ab + (-iy 9 s ab ]S^. (1.3) 

3,a,b 

This results in the presence of an effective staggered field, g s H, upon the application of a uniform field. Such a 
staggered field, which couples directly to the Neel order parameter, is expected to produce an ordered antiferromagnetic 
moment and a gap which scale with field. 

This idea was developed in detail in Ref. [ ||] where it was found that a staggered Dzyaloshinskii-MoriyaQ'El (DM) 
interaction also contributes a roughly equal amount to this effective staggered field. This corresponds to an additional 
term in the Hamiltonian: 

HDM = £(-1)^ ■ (Sj-! x Sj). (1.4) 
3 

It was found that several aspects of the experiments could be explained in detail by this model. These include the 
field orientation dependence of the gap and its scaling with field magnitude. Much of this work used the bosonization 
technique which maps the problem onto the sine-Gordon model, for which various exact results are available. The ex- 
citations observed in neutron scattering were identified with the soliton, antisoliton and "breather" (soliton-antisoliton 
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boundstate) spectruirupf the sine-Gordon model. Additional results, further supporting this approach, were obtained 
by Essler and TsvelikJj The purpose of this paper is to provide more details and some extensions of the results in [ ^| . 
While Cu Benzoate is the only example of such a system that we discuss in the present paper, it should be possible 
to apply our theory to other quasi-one dimensional system with similar crystallographic structure. 

In the absence of a staggered field, the critical behaviour of the antiferromagnet is det erm ined by 3 field dependent 
quantities: the magnetization, m(H) [which determines the soft wave- vectors via Eq. ( [1.21) 1, the spin-wave velocity, 
v(H) and the boson compactification radius (which determines the critical exponents), R(H). All three of these 
quantities can be determined very accurately by numerical solution of Bethe ansatz equations. Furthermore, we 
derive an exact relationship between these three functions using field theory arguments. We also derive the logarithmic 
dependence of R on H , as H — > 0, from the existence of a marginally irrelevant operator, using the renormalization 
group. 

We determine the scaling of gap with field, order parameter with field and susceptibility with temperature. The 
calculations are done including logarithmic corrections which arise from a marginally irrelevant operator and take into 
account both uniform and staggered fields. Furthermore, the exact amplitudes of the scaling functions are determined 
using a recent result of Lukyanov and Zamalodchikov.E3 (After this calculation was finished we received the preprint [ 
0] which gives the same result for the gap, without a discussion of logarithmic corrections.) We also give some further 
discussion of the structure factors S a (q,u>), measured in neutron scattering. We discuss a hidden SU(2) symmetry of 
the model. We prove that the longitudinal structure factor (for "a" corresponding to the uniform field direction) gets 
contributions from only the soliton and anti-soliton intermediate states, in agreement with experiment. On the other 
hand, the transverse structure function gets contributions only from the breathers. Using the approximate SU (2) 
symmetry, we discuss the relative intensity of the various single-particle peaks in the neutron-scattering cross-section, 
taking into account the polarization dependent factors which arise from Fourier transforming the dipole interaction 
between neutrons and spins which were omitted in Ref. [ A comparison is made with experimental results. In 
particular, the problem of determining a consistent value for the DM vector, D is discussed. The susceptibility of 
the sine-Gordon model is calculated, using the integrability of the model, giving essentially the field and temperature 
dependence of the staggered susceptibility. 

In Section II we discuss the DM interaction and the mapping of the system into a Heisenberg model with orthogonal 
uniform and staggered fields. In Section III we treat this problem using conventional spin-wave theory. In Section IV 
we discuss bosonization in the presence of a uniform magnetic field. In Section V we extend the bosonization approach 
to the case with staggered field and analyze the induced gap. In Section VI wc discuss structure factors and compare 
with the observed neutron scattering cross-section. In Section VII we present estimates of the DM interaction based 
on several experimental results. In Section VIII we discuss the magnetization and susceptibility. 



II. EFFECTIVE HAMILTONIAN 

The crystal structure of Cu benzoate is shown in Fig. (|l]) and Fig. (Q). The chain direction is the c-axis. Note that 
each Cu atom is surrounded by 6 ligands with a local symmetry which is almost tetragonal. However the principal 
axes for this tetragonal symmetry alternate along the chain, with the two inequivalent c-axes being rotated by 10° 
relative to each other. These corresponding c-axes are labelled I and II in Fig. (|). The b-axes are the same for both 
Cu sites. Neither of these sets of principal axes correspond to the crystal axes. It is expected that the principal axes 
for the gyromagnetic tensor will also alternate, corresponding to the local tetragonal axes around each Cu ion. The 
principal axis for the anisotropic exchange interaction is expected to be the c' axes, the perpendicular bisector of 
I and II axes. On the other hand, the principal axis for the dipole interaction, which is of roughly the same order 
of magnitude is essentially the c-axis. Combining these two types of contributions to the nearest neighbor spin-spin 
interaction, gives a principal axis which roughly bisects d and c, and is denoted c" in Fig. (^). It is convenient to 
refer the g-tensor to this a"-b-c" coordinate system. From electron spin resonance (ESR) measurements J3 it takes the 
form: 

2.115 ±0.0190 0.0906 \ 
±0.0190 2.059 ±0.0495 = g u ± g s , (2.1) 
0.0906 ±0.0495 2.316 / 

with the ± referring to the 2 inequivalent Cu sites. g u and g s are the uniform and staggered parts of the g-tensor. 
This staggered g-tensor produces an effective staggered field, ±g s H, while the uniform g-tensor produces an effective 
uniform field g u H . In the special cases where the applied field is along the b axis or in the a"c" plane the effective 
staggered field is perpendicular to the applied field and also to the effective uniform field. For general directions of 
the applied field they are almost perpendicular (to within a few %). 
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As discussed by DzyaloshinskiiH and MoriyaS, in magnetic crystals of low symmetry an additional antisymmetric 
exchange interaction occurs, the DM interaction : 

H D M=J2^-(SjxS j+ i). (2.2) 

j 

The possible values of the DM vector Dj can be limited by considering crystal symmetries of Cu benzoate. Firstly, 
the compound is invariant under a translation along the c-axis by two sites. This means the DM vectors are the same 
among the even (or odd) links, but even and odd DM vectors can be different. Secondly, there is a symmetry under 
rotation by angle it about an axis parallel to the b-axis that rpasses through the mid-point of two neighboring sites (J 
and j + 1), along the chains (c-axis). As noticed by Moriya,B this implies the DM vector for the interaction between 
j and j + 1 must be orthogonal to the b-axis. This can be shown as follows: assume we have the DM interaction 
D ■ (Sj x Sj+%). Now apply the rotation described above. It acts on the spin operators as S^' c — > — Sjj+i-fc and 
S%, — > S^j+i-fc - Thus, for the 6-component of D, the DM interaction would be inverted while it is unchanged for the 

a, c-component of D implying that D b = 0. Finally, the crystal structure is invariant under the combined operation 
of one site translation along the chain (c) direction and reflection in the ac-plane. Considering the fact that the spin 
vector Sj is an axial vector, the operation acts as S^'° — > — S^'^ and S| — > <S£ +1 . Since the DM vector is orthogonal 
to the &-axis (and thus one factor of S b always appears in the outer product), the DM interaction term is inverted 
by t he combined operation: D ■ (Sj x Sj+i) — > — D ■ (Sj+i x Sj+%). Thus, the DM vector is alternating as in Eq. 



(1.4). There are apparently no other restrictions that can be placed on the DM vector using symmetry alone. By 
considering a tight-binding model for the exchange interactions it was cstimatccH that D j J is of 0(5g/g) where 5g is 
the deviation of g from twice the identity matrix. 

Apart from the antisymmetric DM interaction, the remaining exchange anisotropy is believed to be quite negligible 
(about 1% of J) and we will henceforth ignore it. Taking D cx z, we may write the Hamiltonian: 



H=\ YpStj-iSv + J*S+S^ J+1 + (h.c.)] 

3 

[S% 7 -1^2 i + SZjSHj+lb (2-3) 
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where J = J + iD. Performing a rotation! of the spins by an angle ±a/2: 

S+ - S+e M/2 , S+ +1 - S+ +1 e- fa /a, (2.4) 

where 

tana = L>/J, (2.5) 
the Hamiltonian is transformed to the standard xxz model: 

H = J2l JS 3 S Ui + ]j Y( S l S I + i + h - c -)l ( 2 - 6 ) 

3 

With some assumptions this anisotropic exchange may cancel the pre-existing one. In any event, it is small and we 
will ignore it. 

Now consider an ext ernal magnetic field, approximating the g-tensor as 2 times the identity matrix. The spin 



redefinition of Eq. (2.4) introduces an effective staggered field. For example, for a uniform field in the x-direction: 

3 3 

Combining the actual form of the g-tensor with the DM interaction, we can obtain effective uniform and staggered 
fields corresponding to an arbitrary applied one. Writing the rotation matrices by ±a/2 around D as: 

U 3 (±a/2) = K u ±K S , (2.8) 

the effective uniform and staggered fields are defined by: 
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H u = [K u g u + TZ s g s ]H 
H s = [n s g u +U u g s )H. 



(2.9) 



In general H s is nearly orthogonal to H u . For small g s and D/J, the staggered field can be approximated as 

H s -9 s H+^jDx {g u H), (2.10) 

namely the sum of two contributions. 

Henceforth, since we are ignoring the small residual exchange anisotropy and assuming that H u _L H s , we will take 
H u to be in the z-direction and H s to be in the x-direction and refer to them as simply H and h respectively. Also 
setting 2/is = 1 we arrive at the simple effective Hamiltonian: 

Heff = J2l J & ■ £+1 - HS t - H-iyS?}, (2-11) 

i 

with h « H . [Note that we have switched the directions of the uniform and applied fields relative to our earlier 
papertj which, unfortunately, contained some inconsistencies of notation.] 



III. SPIN- WAVE THEORY 



In this secti on we summarize the results of spin- wave theory (leading order 1/S expansion) for the effective Hamil- 
tonian of Eq. (2.11). [As far as we know, spin-wave theory results for this problem were first published in Ref. [ 
Although this misses certain features caused by quantum fluctuations in one dimension, it is still quite instructive. 

The classical groundstate is a canted antiferromagnetic structure, shown in Figure (|5|). The spins on both sub- 
lattices lie in the xz plane canted towards the z-axis by an angle 9 from the ± x-axis. The classical energy of this 
state is: 



E{6) /L = -JS 2 cos 29 - hS cos 6 - HS sin 9. 
This is minimized for 9 the solution of: 

4JS 2 sin 9 cos 9 + hS sin 9 -HS cos 9 = 0. 

For h = 0, the solution is: 

sin 9 = H/4JS. 



(3.1) 
(3.2) 
(3.3) 



In order to do a systematic 1/S expansion, it is convenient to regard H and h as being of O(S). We assume that H 
is less than the saturation field, AJS. The leading order spin-wave expansion for a spin pointing in the x-direction is: 



(a] + aj),i\ -(a] - a t ) 



(3.4) 



Here Oj is a boson annihilation operator. To consider small fluctuations about the canted structure we simply write: 



5' 2 



#2i+l ~ TZ-zTZyS^+x- 



where 1Z V is a rotation about the y-axis by — 9: 




(3.5) 



(3.6) 



and 1Z Z is a rotation by ir about the z-axis: 
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n z = 




(3.7) 



Using the facts that: 



n T = n z 



(lZ Z lZy)3i 



-(H y )ii, 



the Hamiltonian may be written in a manifestly translationally invariant way: 

h = ]T[js* ■ n^n z n y s° +1 - H(n y s°) z - h(n y s%] 



(3.8) 



(3.9) 



It is the translational invariance of Eq. (3.9) which motivated the somewhat peculiar looking choice of transformation 
matrices in Eq. ( |3.5| ). Substituting Eq. ( |3.5| ) into the Hamiltonian of Eq. ( |3.9| ), we find that the term of 0(S 2 ) is a 
c-number and the term of 0(S 3 ^ 2 ) vanishes. The term of O(S) is: 



H ^2[(2JScos29 + Hsin9 + hcose)a i j a j + (JS/2)(cos26> - l){a]a j+ i + aj +1 aj) 



+ (JS/2)(cos29 + l)(a f a] +1 + ajOj+i)] 



(3.10) 



Fourier transforming and performing a Bogoliubov transformation that mixes with a)_ k , we obtain a single band 
of spin waves in the paramagnetic Brillouin zone, —tt < k < tt with dispersion relation: 



E(k) = { [2JS cos 29 + Hsm9 + hcos9- JS(1 - cos 29) cos fc] 2 
-[7,5(1 + cos 2(9) cosfc] 2 } 1/2 . 



(3.11) 



The above transformation has allowed us to obtain a single band in the paramagnetic Brillouin zone. We may 
equivalently fold the dispersion relation into the antiferromagnetic Brillouin zone, —tt/2 < k < ir/2. This gives us 2 
branches of spin- waves with dispersion relations: 



E±(k) = {[2JScos26 + Hsin6 + hcose±JS(l~cos26)cosk}< 
-[JS(l + cos2<9) cos k] 2 } 1 ' 2 . 



(3.12) 



While this is a fairly simple and explicit formula for the energies in terms of 9, it must be borne in mind that 9 is 
determined in terms of J, H and h by Eq. (3.2). In the special case h — 0, using Eq. (3.3), we obtain: 



E±(k) = 2JS 

Note that at k = 0, or equivalently k = tt, 



sin 2 k + 2 



— V 

us) 



(cos 2 k ± cos fc) 



1/2 



(3.13) 



E_ = 
E+ = H. 



(3.14) 



The £L mode is the Goldstone mode corresponding to a uniform precession about the z-axis. A non-zero h gives thi s 
mode a gap, pinning the spins along the x-axis. We may calculate this gap, to lowest order in h using, from Eq. (3.2) 
9 = sin" 1 {H/AJS) + 59 where 



Hh 



16J 2 S 2 ~H 2 ' 



To linear order the gap is given by 



A 



dE 2 



59 



dE 2 
~~dh 



■h. 



(3.15) 



(3.16) 
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thus: 



A » y/4JSh[l + (H 2 /8J 2 S 2 )][l - (H/AJS) 2 ] 1 ' 4 + 0(h 3 / 2 ). (3.17) 

Note that A is a singular function of the staggered field, exhibiting a mean field exponent of 1/2. On the other hand, 
it depends only weakly on the uniform field, being almost independent of H up to H of O(JS). While this mean field 
exponent changes when one-dimensional quantum fluctuations are taken into account, it is reasonable to expect this 
weak dependence on H to remain true. The upper mode, E + (0) depends strongly on H but only weakly on h. For 
h « H E+(0) w H. 



IV. BOSONIZATION FOR STAGGERED FIELD 



In the one-dimensional case, an exact picture of the low-energy behavior can be obtained using bosonization and 
RG arguments. Here we summarize the results for the case of a uniform field, but no staggered field. 

We begin with the case where the uniform field also vanishes. The low energy degrees of freedom of the quantum 
spin variables can be represented in terms of a free boson with Lagrangian density: 



£ = \i(M) 2 -v 2 (d x( f ) n 



(4.1) 



[Here v s is the spin- wave velocity which we will generally set equal to 1.] The boson field, <j) can be separated into 
left and right moving terms: 



4>(t, x) = <j> L {t + x) + 4> R (t - x). 



Their difference defines the dual field: 



5R- 



(4.2) 



(4.3) 



For H = there are low energy degrees of freedom at wave- vectors and n. The spin operators can be approximated 
as: 



j ~ 2ttR dx 



+ (-l) J cos^ 



i ■ constant 



J(2ttB.c/>+<P/R) , i(2^R<t>-<t>/R) 



(4.4) 



(The first constant above is universal — The next three (real) constants are not, but C has been recently determined 
using the integrability of the modeO'tj and will be discussed in the next section.) 

For the H=0 Heisenberg model, R = l/y/2n. For the xxz model R (and C) vary with the anisotropy parameter. 
Writing the Hamiltonian: 



H — J~^^[Sj Sj +1 + SjSj +1 + 5SjSj +1 ], 
j 



(4.5) 



exact Bethe ansatz results determine: 



271-ir 



(4.6) 



R varies between l/y/2n and along the xxz critical line, — 1 < S < 1. 

In order to understand the vicinity of the isotropic antiferromagnetic point, <5 = 1, it is convenient to use non-abelian 
bosonization: 



(Jl + Jr) + constant(— l) J tr(ag). 



(4.7) 



Here g is the SU(2) matrix field of the Wess-Zumino-Witten, k=l non-linear a- model (WZW model). Jl and Jr are 
the left and right moving conserved currents associated with the SU(2) symmetry of the spin chain. By comparing 
Eqs. (4.4) and (4.7) we may read off the correspondences between the WZW fields and the free boson fields. These 
are: 



G 



8vr V8tt 



5 « ( _ e iV5F* ) ■ ( 4 - 8 ) 



Here d± = d t ±d x . 

The bosonized spin-chain Lagrangian contains, in addition to the free boson Lagrangian, interaction terms: 



c mt = ^ z J z l Jr + ^{JIJr + «)]• (4-9) 

For the isotropic Heisenberg model, X z = X± = 0(1). Including small anisotropy, X z — X± oc 1 — S. These obey the 
RG equations: 

dX z /d\nE = (4vr/%/3)Ai 

dX^/dhxE = (4vr/V3)A 2 A ± (4.10) 

The RG trajectories are hyperbolas as shown in Fig. (||). For X z > Aj_, they end at the \ z axis, corresponding to 
the xxz critical line. For X z < Aj_ they lead towards strong coupling, corresponding to the easy axis ordered phase. 
Reverting to abelian bosonization, we see that the fixed point Lagrangian contains the extra term: 

^XMJUr - -^A z (O)p t 0) 2 - (d x cj>) 2 }. (4.11) 

Since this is proportional to the free Lagrangian, we can eliminate it by a rescaling of (j>. This corresponds to a 
rescaling of th e pa rameter, R, decreasing it by an amount of O(X z (0)) = 0(\ — 5) in agreement with the Bethe ansatz 
result of Eq. jtq ). 

We now consider the Heisenberg model with a uniform external field (but no staggered field). The extra term in 
the Lagrangian becomes, upon bosonization, 

This term can be eliminated by a redefinition of the boson field: 

<f>(t, x) -> <f>(t, x) + —=x. (4.13) 

V27T 



This leaves the free Lagrangian unchanged. However, it does effect the interaction term and the bosonization formulae. 
The interaction term is changed due to the shift of the ± components of the currents: 



T J- p ±iHx 

J L,R J L,R t: 

j+ j+ p T'iHx 

J L,R. ^ J L,R e 

Jl,R ^ Jl,R (4-14) 

(Note that the phases add, rather than cancel in the interaction term J^Jr-) The effect of these phases on the RG 
equations can be determined from a consideration of the operator product expansion (OPE). One of the OPE's gets 
shifted while the other one does not: 



j£(x)J£(x') 



iH(x'-x) Jl 



J 



j£{x)Jl{x') - — (4-15) 
x — x 

The one-loop RG equations can be conveniently derived using an ultraviolet cut off on the distance between any pairs 
of insertions of the interaction Lagrangian in perturbation theory. (See for example [ pT|.) As the position space 
ultraviolet cut off is increased from a to a', we integrate over that range of separation of the two points, using the 
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OPE. The net effect is that, when the cut-off is smali compared to 1/H the phase factor in the OPE is nearly constant 
and can be ignored. However, when the cut-off is large compared to 1/H the phase factor produces rapid oscillations 
which tend to cancel out t he ter m from the effective renormalization. Reverting to an energy cut-off, this means that 
the RG equations of Eq. ( 4.10 ) are essentially correct for E » H, but for E << H the right hand side should be 
replaced by in the first equation. That is, X z ceases renormalizing, at E of 0(H) whereas Aj_ continues to renormalize 
as before. Thus the RG trajectories are essentially 2 straight lines. For E » H the couplings renormalize along the 
isotropic separatrix, X Z {E) = \±(E), but for E « H \ z is constant and Aj_ renormalizes to 0. See Fig. Thus, 
in order to determine A z (0) we simply need to calculate X(H) using the isotropic RG equations of Eq. ( |4~10| ). For 
H << J, these give: 



X(H) 



V3 



4tt \tl{J/H)' 



(4.16) 



This argument may seem rather naive and it certainly does not give the correct RG trajectory for E of 0(H). However, 
due to the weak, logarithmic, dependence of A on E we expect that this argument gives the correct behaviour of \ Z (Q) 
with H for H « J. This argument then determines the dependence of R on H for H << J: 



2-kB? 



1 



1 



1 



(4.17) 



21n(J/H) 

Precisely this result was obtained from the Bethe ansatzS in the limit H « J. See Fig. @. For larger values of 
H higher order terms in the RG equations would be needed and additional interactions would have to be considered, 
some of which break Lorentz invariance. The net effect is that the H-dependence of R becomes more complicated at 
larger H and also the spin- wave velocity also change with H. We have so far set it equal to 1. It is known to have the 
value it J/2 from the Bethe ansatz, for the Heisenberg model at H = 0. It can be determined numerically from the 
Bethe ansatz integral equations^], and is given in Fig. @. 

Apart from the functions v(H) and R(H) we will also be interested in the behavior of the magnetization, m(H). 
From bosonization we obtain: 



m{H) 



H 

2ttv ' 



as H — > 0. One way of obtaining this result is from calculating the zero field susceptibility: 

2 



X 



K?4 



1 



2ttT 



dx— 

ox 



(4.18) 



(4.19) 



Bosonization leads to an exact relation between the three fun ctions v(H), R(H) and m(H). This follows from 
calculating the susceptibility at arbitrary field, H, using Eq. (4.4) with the corresponding value of R(H): 



dm 
dH 



X(H) 



1 



{2tt) 2 R(H) 2 v{H) ' 

In particular, in the limit of small H this predicts a logarithmic correction to the magnetization: 



H 

2ixv 



1 



1 



2\n(J/H) 



(4.20) 



(4.21) 



Eqs. (4.17), (4.2C) and (4.21) are universal. They should remain true for generic half-integer spin isotropic antiferro- 
magnets in the g apless phase. For the particular case of the nearest neighbour S=l/2 Heisenberg model we can set 
v = 7rJ/2 in Eq. (4.21). Eq. (4.20) agrees very well with our numerical solution of the Bethe ansatz equations. [See 
fi gs. (0), (||), (0).] As far as we know the first (numerical) calculation of these quantities for < H < 2 J were Ref. 



m for m(H), Refs. [ 18, 19 1 for R(H) and Ref. [20] for v{H). (Numerical calculations in Refs. [ pi M were based 



on the Bethe ansatz integral equations of Ref. [ |16 [ 

Once m, v and R are deter mined from the Bethe ansatz, all low energy properties of the system are determined by 
bosonization. From Eq. (4.4) the magnetization can be written: 



m(H) = 



1 f , d± 
2irR(H)J dx' 



(4.22) 
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Therefore the exact formula for the field-induced shift in <f> is: 



2%R(H)m(H)x. 



(4.23) 



All low energy Green's functions are then determined from Eq. (4.4) after shifting 4>{x). In particular we see that G z 
has the soft wave- vectors and ir ± 2ivm(H) whereas G has the soft wave- vectors ±2irm(H) and it. 



V. BOSONIZATION FOR NON-ZERO STAGGERED FIELD 



Now we consider the effective Hamiltonian of Eq. ( [2.11 ) with both uniform and staggered fields non-zero and 
h « H. We begin by using the results of the previous section to obtain the h = theory with the shifted and 
rescaled boson field characterized by R(H), m(H) and v( H). For h = 0, upon making these tranformati ons , the 
Lagrangian density is simply the free boson one of Eq. fl4.1| ). From the bosonization formulae of Eq. (4.4) the 
staggered field adds the interaction term: 



C-int = hC cos (2nR(f>). 
Noting the duality transformation between 4> and (j>: 

d t 4> = d x <p d x 4> = -d t (, 
we may also write the free Lagrangian in terms of <p: 

c = \[(d t 4>) 2 -(dj) 2 }. 



(5.1) 



(5.2) 



(5.3) 



Hence we have the standard sine-Gordon field theory. An impressive array of conjectured exact results are available 
on this model, which can be brought to bear on the spin-chain problem. The interaction term is sometimes written 
as cos f3cj), so we see that we have (3 = 2irR. In the next section we will discuss details of the excitation spectrum of 
this model. In this section we discuss the dependence of the gap on the uniform and staggered field. 

The renormalization group scaling dimension of this operator is irR 2 from which it follows that the gap scales as: 



A 
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A 



J 



l/{2-xR.(H) 2 ) 



(5.4) 



for some function A. Since, for small uniform fields, irR 2 1/2, it follows that the exponent is approximately 2/3, as 



-> at fixed H. This is a reasonable order of limits for 
The exponent is determined by R(H), given in Figure 



found in the experiment. Note that this formula is valid for h 
describing the experiments since h is only a few per cent of H. 

(§■ 

We wish to improve on this result in two ways. First of all, if we consider the case where H is strictly 0, then this 
formula is modified to: 



A/J 



>"'7 



2/3 



hi 



1/6 



(5.5) 



We will calculate exactly the amplitude Aq . This is not of direct relevance to the experiments however since they are 
in the opposite limit H » h. More importantly, we can determine the amplitude function A{H/ J) in Eq. (5.4) for 
H « J: 



A 



Aln 



1/6 



(5.6) 



We will also determine ex actly the numerical factor, A, which is different than Aq in Eq. (5.5). While the logarithmic 
factors in Eqs. (|]^) and (5J3) are universal and follow from an RG treatment of the marginal interaction, the exact 
numerical coefficients are specific to the ordinary nearest neighbour Heisenberg S=li2 model and are obtained by 
using a remarkable exact conjecture made recently by Lukyanov and ZamalodchikovO extended to the Heisenberg 
point following the method in Ref. [ [l|. [See also Ref. [ Il4|.] 

Our calculations follow the notation of Refs. [O] and []13|. We consider the renormalization group equation obeyed 



by the effective coupling constant g{E) multiplying the cos interaction in Eq. (5.1), with ba re v alue hC/J for a bare 
cut off J, taking into acount, to linear order, the effect of the marginal interactions of Eq. (fi~9|). This is: 
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dg 



dQjxE) 

The anomalous dimension is given, to low order, by: 

7 



-[2-7(0>A)]fl. 



1 —A 

2 V3 2 



(5.7) 



(5.8) 



A 2 obeys the RG equation of Eq. ( 4 . 1 0| ) . By the usual scaling arguments, we determined the gap, A by reducing 
the ultraviolet cut-off down to a scale A such that g(A) is 0(1). Taking into account the dependence of the effective 
coupling constant, g(E) on the bare coupling constant, h/ J, then determines the dependence of the gap on h. In the 
case H = 0, the RG flow of X Z (E) is given by Eq. ( 4. 10| ) for all E. On the other hand, for finite H, X Z (E) essentially 
stops renormalizing at a scale of order H. Integrating Eq. ( |5.7| ), gives: 



9(E) = 



-3/2 



-tt/a/3 d In E'\z (E') 



(5.9) 



This integral can be conveniently evaluated by changing integration variables to X±(E') using the second of Eqs. 
(P~10|). This gives: 



g(E) « g{J){E/J)-^[\^E)/\^{J)Y/\ 
Determining the gap by the condition g{A) = 1 and setting g( J) cx h/ J, gives: 

h/J = B(A/J) a ' 2 [4x\ ± (A)/VE\ 1/4 , 

for some non-universal constant B, of 0(1). 

In the case H = 0, the solution of the RG equation Eq. ( 4.10| ) is: 



(5.10) 



(5.11) 



A X (A) 



V3 



Thus, 



4?rln(J/A) 

A/J = B~ 2/3 (h/J) 2/3 [ln{J/h)} 1/6 [l + 0{l/\n{J/h))}. 



(5.12) 



(5.13) 



For non-zero the first of Eq. ( 4.10 ) is only valid for E » H. At lower energies, X z stops renormalizing. Its 
fixed value at low E determines R(H): 



27rA 2 (0)/V3 = 1 - 2ttR(H) 2 



1 



2\n(J/H) 



(5.14) 



We can extend somewhat the accuracy of our results to larger H , by expressing the subsequent results in terms of 
R(H), determined numerically from the Bethe ansatz, rather than by using the above asymptot ic sm all H result for 
R(H). For A << H, we use the second of Eq. ( 4.10 ) with X z fixed at A z (0) as given by Eq. ( 5.14 ) and the initial 
condition 



to obtain: 



X±(H) « X Z (H) » A,(0), 



X ± (E) » —[1 - 2nR(H) 2 }(-) 2 ^ R ^ 2 l 
2tt H 



Setting E — A and substituting into Eq. (5.11 



gives: 



(\-2ttR 2 )/2 



[2(1 - 2-kR 2 )} 



1/4 



(5.15) 



(5.16) 



(5.17) 



Thus we obtain Eq. (5.4) with 
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A(H/J) m { B 



(l-2?ri? 2 )/2 



-1/(2-tt_R 2 ) 



[2(1 - 27Ti? 2 



I 1/4 



(5.18) 



Note that in Eq. (5.4) and (|5.18| ) R is a function of H, sh own in Fig. (|J) and given approximately by Eq. ( 5.14 ). As 
H — > 0, we may evaluate -4(i/) explicitly from Eq. ( |5.14 ). Using: 



l/41n(J/ff) 



= 1/4 



this gives: 



A{H) -> B- 2/3 e- 1/6 {\n(J/H)) 1/e . 



(5.19) 



(5.20) 



Note that the same numerical constant, B, occurs in the H = case, Eq. ( 5.13 ) and the H >> A case, Eq. (5.4) 
and (5.20). However, in the latter case it gets multiplied by an extra factor of e -1 / 6 



Finally we wish to determine the dimensionless amplitude, B, appearing in Eq. (5.11) and below. This can be done 
using two remarkable results. One of them is the ex act proportionality constant in the bosonization formula for the 
staggered part of SJ, that is the constant C in Eq. (4.4). This fixes the coupling constant in the sine-Gordon model, 
Ch, in Eq. (|]l]). The other recent result is the exact relationship between the sine-Gordon coupling constant and the 
mass of the lightest particle in the spectrum of the sine-Gordon model, which is the gap, A. This determines the exact 
relation ship between A and h. This calculation was done in zero uniform field for the xxz S=l/2 antiferromagnet 
of Eq. (4.5). The calculation was performed for all 5 along the critical line, — 1 < 5 < 1. Due to the logarithmic 
corrections at the isotropic point, 5 — 1, an ad dition al calculation is needed at that point. This can be done using 
the RG. We essentially just need to apply Eq. (5.11) to the case of H = but S slightly less than 1. Comparing to 
the exact result for all 6 < 1 then determines the coefficient, B. 

For S < 1 and H = we use the RG equations of Eq. ( 4.10p at low energies. The RG flows are hyperbolas 
terminating on the positive X z axis. These flows are conveniently labelled by: 



cos 1 S 



e = 4tt\ z (0)/V3 = 2[1 - 2ttR 2 } = 2 

The solution of Eq. ( 4.10| ) with this initial condition is: 

4ttA_l(A) _ e 

V3 ~ sinh[eln(J/A)]" 



Substituting into Eq. (5.11) gives: 



J \jj lsinh[eln(J/A)] 



1/4 



(5.21) 



(5.22) 



(5.23) 



Note that in the isotropic limit, e — * 0, we recover our previous logarithmic result of Eq. ([5.5|). On the other hand, 
taking A/ J — > with e held fixed it gives: 



h /A x3/2+e/4 



(5.24) 



By comparing this to the exact result for arbitrary <5 (and hence e) we may extract the value of the amplitude, B. 
In Ref. [ 12 1 the spin correlation function in the xxz antiferromagnet is shown to have the asymptotic behavior: 



< S$S? >^ (-1) 



■C(Rf 



with an exact expression determined for the amplitude C(R) 2 . In the isotropic limit, R — > 

C{R) 2 1 
2 4e 1 /2 7r 3/2 ■ 



(5.25) 



(5.26) 



In Ref. [ 12| the operator cos2ttR(/) is normalized so: 
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(cos[27ri?0(O)]cos[2vri?0(r)]) -► -\r\- 2wB? 



(5.27) 



(after accounting for a difference in normalization of the free boson Lagrangian by a factor of 87r) . This determines 
the exact proportionality constant in the bosonization formula of Eq. (4.4): 



S* sa (-l) j C(R) cos[2ttR$\. 



(5.28) 



Hence the coupling constant in the sine-Gordon Lagrangian is precisely C(R)h. The exact relationship between this 
coupling constant and the mass, A, of the soliton of the sine-Gordon model is: 



C(R)h 



.2-nR" 



T{ttR 2 /2) 



v 1-kB? ^(1-7^2/2) 





v 2(l-7rflV2)) 




2r 


( ttR 2 \ 




\4-2ttB: 2 ) 



2—nR A 



(5.29) 



Here we have inserted, by dimensional analysis, the spin-wave velocity v. Taking the isotropic limit on both sides of 
this equation, using v — ttJ/2, gives: 



h 



A 



2V2 7 r 3 / 4 e 1 /4 j V J 



3/2+e/4 



r(i/4) 



27r 3 / 4 r(3/4) 



£(2/3) 

r(i/6) 



-i 3/2 



We see that this is equivalent to the RG result of Eq. (5.24) with the amplitude determined to be: 



B = 2 J/ 



i /4 r(i/4) 



r(3/4) 



r(2/3) 



r(i/6) 



3/2 



From Eq. (5.13), for H = 0, the gap behaves as: 

A/J-> A)(V^) 2/3 (ln(J/^)) 1/6 , 

with: 



A = B 



-2/3 



,-1/6 



r(3/4) 



r(i/4) 



2/3 



£(1/6) 

r(2/3) 



1.77695. 



(5.30) 



(5.31) 



(5.32) 



(5.33) 



Since we were not aware of the result [ Q at the time, in [ || the behavior of the gap with staggered field (for H = ) 
was estimated numerically by extrapolating Lanczos results for lengths up to 22 sites. A very good fit to Eq. fl5.32j ) 
was obtained with: 



A) ~ 1-85. 



(5.34) 



Considering the nu merica l difficulties related to logarithmic corrections this is remarkably good agree ment with the 
exact result of Eq. ( 5.33 ). In the experim enta lly relevant case, A < < H, the gap behaves as in Eq. (5.4) and, for 
small H/J, the amplitude is given by Eq. (5.6) with, from Eq. (5.2C), 



A = B- 2/3 e- 1/6 w 1.50416. 

Thus our expression for the gap becomes: 

A/ J -» 1.50416[ln(J/i?)] 1 / 6 (VJ) 1/[2_7rfl(ff)2)1 - 
For larger H/J greater accuracy might be obtained by using Eq. ( |j.1S ) with B given in Eq. ( 5.31 ). That is: 

-1/[2-kR(H) 2 ] 



A/ J -» |.422169(J/i?) [1 - 27r ' R(ff)2l/2 [2(l - 27ri?(i?) 2 )] 1 / 4 }~ 



(h/J) 



l/[2-nR(Hy] 



where R(H) is given, from the Bethe ansatz, in Fig. (g). Inserting its asymptotic value at low H, 

1 



2nR(HY 



1 



(5.35) 
(5.36) 

(5.37) 
(5.38) 



2\n{J/HY 

gives back Eq. ( 5.36| ). We note that, to actually fit the experimental data, we take h = cH for some constant of 
proportionality which depends on field direction but is generally of order a few %. Thus the actual scaling of gap 
with field is not a pure power law. 
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VI. STRUCTURE FUNCTIONS 



For H = it is convenient to use non-abelian bosonization, Eq. (4.7) so that the interaction term is written: 

£ mt oc tr{ga x ). (6.1) 

In this case the model has an SU(2) symmetry. Note that when both uniform and staggered fields vanish and ignoring 
the marginal operator, the symmetry is actually SU(2)xSU(2): 



9 



UgVK (6.2) 



These two independent SU(2)'s act on left and right-movers separately. The ordinary SU(2) symmetry of the spin 
chain is the diagonal subgroup with U = V. This symmetry is broken by the staggered field. However, a different 
SU(2) subgroup of the original SU(2)xSU(2) survives for which: 

V = a x Ua x . (6.3) 

We may redefine the field g by 

9 - go' (6.4) 

in which case the interaction becomes trg, which has the diagonal SU(2) symmetry. In fact, this continuum limit 
interaction arises from a staggered Heisenberg exchange interaction, as occurs in the spin-Peierls problem. The 
equivalence of the continuum limit of these two apparently very different problems, is a non-trivial consequence of the 
chiral symmetry which maps (—l) J Sj into (— iySj ■ Sj+i. 

If we assume that the only effect of the uniform field is to add a term to the Hamiltonian: 



then the uniform field can be removed by the gauge transformation: 



<f>(x) -> 4>(x) + —f=x, (6.6) 

V 27T 



or equivalently: 



JI,r^Jl,r-^- (6-7) 



[Here we have set v = 1.] This transforms the matrix field g as: 

g(x) -> e lHxa */ 2 g(x)e lHxa * /2 . (6.8) 

This gauge transformation leaves invariant the staggered field term trga x . Thus the exact SU(2) symmetry remains 
in this approximation. However, additional irrelevant terms in the Hamiltonian in the presence of a uniform field 
break the exact SU(2) symmetry, as evidenced by the change in the parameter R with field. The SU(2) symmetry is 
only present for R = l/v2~7r. 

The non-abelian bosonization formula for the staggered part of the spin operators (for H = 0) is: 

Sj m (-lyCtrga. (6.9) 



The magnetic field leads to the gauge transformation of Eq. (6.8). The parts of the spin operators with wave- vectors 
near tt thus become: 

Sj w Ccos(7T7')trg(T a (a = x, y) 

S3 C{e l ^ +H ^trg{\ + ct 2 )/2 - e i( -^ H ^trg(l - a z )/2} (6.10) 

Now making the transformation g — ► ga x , this becomes: 
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« Ccos(7rj)trg 
Sy m iCcos(7rj)tr3cr 2 

S| m C{e l ( 7r+ ^'tr ff a-/2 - e l{7 ^ H) Hxgn+ /2}. (6.11) 

In general the spectrum of the sine-Gordon theory consists of the soliton, anti-soliton and breathers (soliton-anti- 
soliton boundstates)£ll In the SU(2) symmetric case, the excitation spectrum of the sine-Gordon model with f3 2 = 2tt 
consists of a triplet composed of soliton, anti-soliton and lowest breather and a second breather, heavier by a factor of 
V3- The degeneracy of the triplet is a result of the SU(2) symmetry. Due to the SU(2) symmetry, the 3 elements of 
the triplet are produced by the 3 operators trga with equal intensity, whereas the singlet is produced by the operator 
teg. Thus setting H — (but not h) the structure functions G vv and G zz would be equal. Note that, for H = 0, 
S z oc tvga v , creates the y-polarized member of the triplet. This can be regarded as a linear combination of the soliton 
(created by tiger - ) and the anti-soliton (created by tiga + ). G zz (it,uj) consists of 2 identical contributions from the 



a~ and <j + terms in Eq. (6.11). Each contributes exactly (1/2)G vv (tt, u>). The effect of H is to split the soliton and 
anti-soliton contributions to G zz into two separate contributions at different wave- vectors ir ± H . Thus, ignoring the 
small SU (2) symmetry breaking: 

G zz {-k±H,uj) = (l/2)G yv {n,uj). (6.12) 

It is also interesting to note that the staggered part of the energy density is given by: 

Sj ■ S J+1 oc (-lytrg. (6.13) 

This operator couples to lattice displacements (phonons) and is used to describe Raman scattering experiments. Upon 
making the gauge transformation of Eq. (|6.8[) and the redefinition of Eq. (|6.4|) this becomes: 



S 3 ■ S j+1 oc e *C*+*)itr^- + e *(^b tr ^± (6 . 14) 

Thus this operator also creates the soliton and anti-soliton. Hence this theory predicts a single particle excitation 
observable in Raman scattering at the same field-dependent wave- vector and frequency as the incommensurate mode 
observed in neutron scattering. 

Upon allowing for SU(2) symmetry breaking the radius changes. After making the gauge transformation, a U(l) 
symmetry still survives, corresponding to shifting <p by a constant. The triplet is now split, with the lowest breather 
having a different mass than the degenerate soliton antisoliton pair. Since the operators e ±l ^^ R have charge ±1 with 
respect to this U(\) symmetry, we see that the soliton and anti-soliton are created by the q = 7r± 27rm Fourier modes 
of Sj respectively. The breathers can be classified as even or odd with respect to the discrete symmetry <p — > — <f>. 
The odd breathers are created by the q = it component of S v and the even breathers by the q — 7r component of 
S x . It can be shown that even and odd breather alternate in the spectrum of the sine-Gordon model. Furthermore, 
the number of breathers increases with decreasing R. A third breather drops below the soliton antisoliton (ss) 
continuum immediately as soon as R decreases below the isotropic value, 1/V2tt with another one dropping below 
the ss continuum each time 2/ttR 2 passes through an integer. The mass of the n th breather, expressed in terms of 
the soliton mass, M, is: 

M n = 2Msin(n</2), (6.15) 

where: 

1 2 



£ ttR 2 



1. (6.16) 



Thus the odd-numbered breathers contribute single-particle poles to G vy and the even-numbered ones to G xx while the 
soliton and antisoliton contribute single-particle poles to G zz . In addition various multi-particle continua contribute 
to the three spectral functions. 

For a field of 7 T. we estimate ttR 2 — .41. There are 3 breathers at this point, with masses .79M, 1.45M and 
1.87M. A resolution limited peak was observed at q = 1.227T, at energy .22 meV. We identify this with the soliton (or 
anti-soliton) contribution to G zz ; hence M=.22meV. A resolution limited peak is clearly observable in the neutron 
scattering data at q — tt and an energy of .17meV=.77M. This agrees very well with the prediction for the first 
breather mass. 
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We may also test the SU(2) prediction of Eq. (6.12). The SU(2) symmetry is broken by various small effects as 
exemplified by the fact that ttR 2 ^ 1/2. In particular, this implies that G vy (ir,u>) has a second peak, of very low 
intensity, corresponding to the third breather. Ignoring these effects we expect the intensity of the lowest breather 
peak in G yv to be approximately twice the intensity of the soliton peak in G zz . 

However, before a comparison can be made with experiment it must be taken into account that the unpolarized 
neutron scattering cross-section contains an important direction dependence arising from the Fourier transform of the 
dipole-dipole interaction between the neutron and the spins. The cross-section can be written: 

a(k, u) = ]T(1 - k 2 a )G aa (k, u)f{k), (6.17) 

s 

where k is a unit vector in the direction of k and the function f(k) is slowly varying. Thus the soliton, even breathers 
and odd breathers are weighted by different factors 1 — k 2 , 1 — k 2 and 1 — k 2 , respectively. We also have to consider 

the variation of f(k) in examining the relative intensity of solitons to breathers since they occur at different values of 
k z . It must be recalled that x here refers to the direction of the effective staggered magnetic field (and z the direction 
of the uniform field). In the neutron scattering experiments the field was along the b-axis. We note that the a" axis 
is rotated by about —17° from the crystal axis, a. [We define the rotation angle in the ac- (a"c"-) plane so that c 
axis is +90° rotated from a-axis.] 

Strictly speaking, we must take into account the effect of the redefinition of the spin operators, discussed in Section 
II, that was used to eliminate the DM interaction. Letting Sj to be the rotated spin operators, defined in Eq. (p~|). 
and inverting the transformation, we may write the structure function for the original spin operators in terms of the 
structure function for the rotated operators, which we write as G ab (k). In this way we obtain: 

G xx {k) = cos 2 (a/2)G xx (k) + sm 2 (a/2)G vv (k - tt) 

G vv (k) = cos 2 (a/2)G yy (k) + sm 2 (a/2)G xx (k - tt). (6.18) 
G xy remains due to translation invariance and the G az are unaffected by the transformation. Here x and y refer to 



two axes orthogonal to D ( not orthogonal to H as in most of this paper.) We expect the second terms in Eqs. (6. IS 



to be negligible since a is small. Furthermore, the G aa {k) are small for k « also making the second terms in Eq. 



( 3.18 ) smaller than the first for k near tt. Henceforth we ignore this small correction and simply use G aa (k) « G aa {k). 

The direction (and magnitude) of the effective staggered field depends on both the staggered part of the ^-tensor 
and DM interaction. Since the DM interaction in Cu Benzoate is unknown, the direction x is not known at present. 
(The DM interaction in Cu Benzoate may be estimated from various experimental results based on the present theory, 
We will discuss this issue in the next section.) However, the direction x can be deduced from the polarization analysis^ 
of the neutron scattering experiment. Dender et al. analysed the polarization of the neutron scattering at constant 
energy Tilo = 0.21meV and various momentum transfers. For magnetic field H = 7T || b and momentum transfer 
along the chain tt, this should probe the lowest (n — 1) breather. As discussed above, this odd breather is polarized 
orthogonal to the total effective staggered field. In Ref. [ |l| it was claimed that the observed scattering is polarized 
in the a" direction. Also note that a misstatement of the crystal orientation occurred in Ref. [ |) so that the wrong 
sign appeared there for k ■ a.@ Correcting this error, the polarization is +18° from a-axis, or equivalently +35° from 
a"-axis. This implies that, for an applied field in b-direction, the total effective staggered field —72° (+108°) from 
a-axis, or equivalently —55° (+125°) from a"-axis. [There is actually another feature of the polarization analysis in 
Ref. [ |^] which appears inconsistent with the theory presented here, namely the polarization analysis in zero field. 
The strong dependence of the intensity on the a-component of the momentum is taken to indicate that G bb « 0. Our 
work ignores any anisotropy in the zero field limit and therefore predicts G bb — G aa — G cc in that limit. We do not 
understand the source of this discrepancy at present.] 

The constant-Q (momentum transfer) scan experiments, sensitive to the breather modes, were carried out at the 
fixed momentum transfer: 

(k-a,k-b,k- c)/2n = (-.3, 0, 1), (6.19) 

where the lattice constants are a = 6.91A, b = 34.12A and c = 89. 3A Here we have again corrected the misstatement^ 
of the crystal orientation in Ref. [ mentioned above. Note that the antiferromagnetic wave-vector, in the chain 
direction, is actually 27r/c rather than the normal tt/c, because there are 2 Cu atoms per unit cell along the c-axis. 
The crystal axes a, b and c are essentially orthogonal. Thus, in the a-b-c system: 

k = (-.26,0,-97), (6.20) 
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k is rotated +105° from the a axis, or +122° from the a"-axis. We note that, this k is almost parallel to the direction 
of the total effective staggered field estimated from the polarization analysis above. 

The constant-Q scan for the soliton modes was done for a slightly different momentum transfer 

(k-a,k-b,k-c)/2TT= (-.3, 0, 1.12). (6.21) 

However, the direction of this momentum transfer k is almost the same as the above and we will ignore the difference. 

The fact that, an intense first breather peak is observed experimentally, supports the deduction that k is nearly 
parallel to the total effective staggered field. In fact, assuming that k is completely parallel to the staggered field, 
the approximate SU(2) prediction is that the lowest peak at k = tt (from the first breather in G yy ) should have twice 
the intensity of the lowest peak at k — tt ± H (from the soliton/antisoliton in G zz ). This prediction is only a very 
approximate one due to the breaking of SU(2). There could also be corrections from the function f(k) in Eq. ( 3.17[ ). 



Experimentally this ratio appears to be about 2.8. This is perhaps satisfactory agreement. This agreement is only 
worsened if k deviates from the direction of the staggered field. 

Thus, analyses of polarization and of the scattering intensity of the lowest breather mode are consistent, and 
appar ently lead to the conclusion that the total staggered field for H||b is almost parallel to the k direction of Eq. 
( |6.20 ) used in the constant-Q scan. 



On the other hand, in Ref. [ |3J we discussed a feature in the experimental data at k — tt, and lu = .34meV — 1.55M. 
This is very close to the predicted mass of the second breather, which contributes to G xx . A recent calculation, based 
on integrability of the sine-Gordon model, indicates that the relative intensity of the second breather in peak in G xx 
should be roughly 1/2 of the intensity of the first breather in G yy . However, we expect this to have essentially zero 
intensity in the neutron scattering cross-section due to the factor of 1 — fc^ in Eq. (6.17). [In Ref. [ [| an apparently 



good agreement between theory and experiment was obtained because the factors of (1 — k 2 ) were not included. 
In Ref. [ |j| intensities were not considered.] A possible resolution of this disagreement is that this "feature" at 
lu = .SAmeV, discussed by two groups of theorists, is just noise. As stated in Ref. [ ^] "Given the quality of the data, 
this double-gap conjecture is highly speculative" . Clearly more data is needed to determine whether or not there is 
really another sharp peak at this frequency. An even more statistically insignificant feature in the data, at k = tt and 
lu = AAmeV, was discussed in Ref. [ ^| where it was interpreted as the third breather peak in G yv . Both the energy 
(« 2M) and the intensity (very approximately 1/6 the intensity of the first breather) agree with the predictions of 
the sine-Gordon model. Note that in this case the factor of 1 — kz is common to first and third breathers so it doesn't 
affect the intensity ratio. However, once again, considerably more data is needed to determine if there is really a peak 
at this frequency. 

Another striking feature of the neutron scattering data is a second resolution limited peak at k = tt and lu — 
.8meV ~ H. This also has a possible interpretation in our field theory approach. It is natural to assume that this 



peak actually comes from G zz . At this wave-vector, from Eq. (6.11), this is proportional to < tvg<j^tvga + > at 
wave-vector k = H. We expect the continuum limit to hold for some range of wave-vectors close to the gap minimum 
at wave- vector tt + H . Thus, at least for weak enough fields, we would expect the soliton to persist as a single-particle 
excitation in G zz up to wave-vector tt. Its energy should obey the Lorentz invariant formula: 

u) = \/M 2 + k 2 = \JM 2 +H 2 . (6.22) 

(Here k is measured from the incommensurate wave- vector k + H and v is set equal to 1.) Since H » M this gives 
approximately u> = H, as seen in the experiment. The intensity of this feature in G zz can be easily calculated. The 
result follows from the fact that the k near tt parts of the spin operators are all Lorentz scalars. The matrix element 
between groundstate and a single particle excited state of a Lorentz scalar operator is independent of the momentum 
of the particle, by Lorentz invariance (assuming a Lorentz- invariant normalization of the state). It then follows that 
the soliton and anti-soliton peaks in G zz have an intensity that is proportional to 1/lu. The energy is approximately 
four times higher at k — tt. We must also take into account that both soliton and anti-soliton are contributing at 
k = tt which increases the intensity by a factor of 2. Thus, we expect the single-particle peak at lu ~ H, k = tt to 
have an intensity approximately 1/2 that of the peak at k = tt + H. Experimentally this ratio looks somewhat larger 
than 1/2 but it must be remembered that the peak at lu ss H, k — tt is sitting on top of a background from G yy . 
Taking this into account, the agreement looks fair. 

There is actually a possible objection to this argument. The same reasoning would seem to imply sharp peaks 
near lj = H at k = tt + H coming from G yy and G xx . These were not observed experimentally; at most a small 
shoulder, was observed beginning at lu = H for k = tt + H . This may simply mean that the breathers have merged 
into the continuum by this wave-vector, (due to non-relativistic effects not contained in the continuum limit field 
theory) whereas the soliton has not. 
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We note that spin-wave theory fails to capture the one-dimensional quantum fluctuation dominated physics in 
various ways, ft predicts a single low energy mode with A oc /i 1 / 2 , instead of h 2 ^ 3 with soft wave-vectors 7r and 0, 
missing the incommensurate shift. It also predicts another single particle mode at energy approximately H, at the 
same wave- vectors. 



VII. ESTIMATE OF DM VECTOR 



In the present framework, the only unknown parameters of Cu Benzoate are the DM vector, which has not been 
determined directly in previous studies. Based on the present theory, we can in principle determine the DM vector 
from several experimental results. Actually, there seem to be no solution that can perfectly fit all the available 
experimental data, as explained below. Presumably, a precise error estimate on an experiment gives a permissible 
region for the DM vector, and such constraints from several experiments would give a region of possible values of the 
true DM vector. However, it should be noted that we have been ignoring the interchain effects, irrelevant operators 
etc. which might be necessary in such a precise discussion. 

Firstly, as argued in Sec. II D must lie in the a" — c" plane, leaving two free parameters. The total staggered field 
is determined by eq. 

Here we list the constraints on the DM vector from various experiments. 



A. Angular dependence of the gap 

As observed in Ref. [ the induced gap is strongly dependent on the direction of the applied uniform field. In the 
present theory, this is accounted by the field direction dependent constants of proportionality between the uniform 
and staggered magnetic fields. As discussed in Sec. II, the proportionality constant should be given by the g-tensors 
and the DM vector D. Since the g-tensors were obtained previouslyLJ, the measured gap can be used to determine D. 

The gap is proportional to essentially the 2/3 power of this effective staggered field. The specific heat measurements 
of Ref. [| were fit by the authors to the specific heat of a collection of free massive relativistic bosons. The masses 
were found to scale approximately as H 2 / 3 with a direction dependent amplitude in the a" : b : c" = .55 : 1.0 : 2.0. 
In Ref. [H), this-ratio was used to estimate the DM vector. 

Very recentlyu the specific heat of the sine-Gordon model was calculated from the thermal Bethe ansatz and fit to 
the Cu benzoate data. Again a good fit of the gap to H 2 ! 3 was obtained for fields in the b or c" directions with a 
somewhat different amplitude ratio b : c" = 1.0 : 2.2. (The velocity, v(H) is another parameter in the fit. This may 
also be determined from Bethe ansatz for the S=l/2 chain in a uniform field. A slightly better fit to the specific heat 
data was obtained in [ Q by letting v(H) be a free parameter.) 

A reasonable fit was not obtained for the field in the a" direction where the specific heat data is nearly linear. This 
suggests that, for H|| a" the apparent gap structure was either due to some sort of experimental error or due to other 
mechanismsQ than the effective staggered field. In any case, it seems that the effective staggered field for H|| a" is 
rather close to zero. This implies the cancellation of the effective staggered field coming from the staggered (/-tensor 
and the DM interaction. It is not quite unnatural, because for the applied field in ac-plane, both the staggered field 
generated by the staggered g-tensor and the DM interaction point to b-direction. Thus there is a direction in the 
ac-plane where the cancellation occurs, for a wide range of parameters. Actually, the cancellation of the staggered 
field is also consistent with the Electron Spin Resonance (ESR) result, which will be explained later. 

Assuming the cancellation of the staggered field for H|| a", and that observed gap for H||b and H|| c" are entirely 
due to the staggered field, the ratio of the proportionality constants between staggered and uniform fields for H|| a", 



H||b and H|j c" are 0:1: (2.2) 3 / 2 = 3.26. The ratio gives a constraint on the DM vector through (2.9) 



B. Magnitude of the gap 

Based on several exact results on the sine-Gordon field theory and on the S — 1/2 Heisenberg antiferromagnetic 
chain, we have determined the magnitude of the gap for a given staggered field h. Thus, comparing this with the gap 
estimated from the specific heat measurement, we can fix the proportionality constant c between the staggered field 
h and the applied field H (h = cH). 

Here we use only the result for H|| c", which is presumably most reliable. For H|| c", the gap is very well fit by 
the power law A = kH 2 ^ 3 , without introducing the logarithmic correction. The proportionality constant is given by 



k = 1.316 if A and H arc measured in units of Kelvin and Tesla, respectively. By comparing this with eqs. (5.32) and 
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( |5.33 ), where we assume the logarithmic factor to be clos e to unity for the present case, we obtain c = 0.111. This 
also gives a constraint on the DM vector through eq. ( |2.9|) . 



C. ESR linewidth 

An anomalous broadening of ESR, which is strongly dependent on the field direction, at low temperatures was 
observcdE3 in Cu Bcnzoatc. The mechanism of this broadening was left unexplained- However, we have recently 
developed a field-theory approach to ESR on quantum spin chains at low temperature.cJ According to the theory, the 
contribution of the staggered field to the ESR linewidth is given by 

roc^. (7.i) 

This diverges at lower temperature, in agreement with the experiment. The direction dependence can also be explained 
by the direction-dependent proportionality constant between the effective staggered field and applied uniform field. 
Actually, this is consistent with the previous discussion on the field-induced gap, at least qualitatively. In particular, 
for H|| a", the low-temperature anomalous part of the ESR linewidth vanishes. This implies the cancellation of the 
staggered field for H|| a". While this appears to contradict the apparent gap found in Ref. [ Q, it is rather consistent 
with more refined analysis discussed in the last subsection. In addition, H|| c" gives the largest linewidth, which is 
consistent with the larger field-induced gap for H|| c". 

On the other hand, the ratio of the staggered field is not quantitatively consistent with the specific heat measurement. 
The estimate of the staggered field is somewhat subtle because there are also contributions to the ESR linewidth from 
other sources (most importantly exchange anisotropy/dipolar interaction). The low-temperature anomalous part, 
which is related to staggered field, appears to be approximately 1 : 4.6 for H||b and H|| c". This gives the ratio of the 
staggered field 1 : 2.1 for H|jb and H|| c". This is smaller than expected from the specific heat analysis. 

D. Neutron Scattering 

As we have discussed in Section VI, the analyses on the polarization and intensity of the first breather suggests 
that the total effective staggered points to —72° (+108°) from a-axis, or equivalently —55° (+125°) from a"-axis, if 
the external field is applied in b-direction. This gives another constraint on the DM vector. 

E. Summary of the estimate of the DM vector 

There are several experimental data which give some constraints on the DM vector, and they are not perfectly 
consistent. The estimate of DM vector is also sensitive to the assumed form of the g-tensor, extracted from ESR 
measurements.l3 The experimental data are presumably subject to several errors which has not been identified precisely. 
We hope that more experimental data will be available in the future to make more precise comparison with the theory. 

The ESR linewidth was measured for various directions of the applied field. Thus it is perhaps quite reliable that 
the staggered field precisely cancels at H|| a" . As we have discussed, this is rather consistent with the refined specific 
heat result. This gives a single constraint on the DM vector. In the linearized approximation, it reads 

h a n = 0.0190 + 0.0453^p - 1.058^p = (7.2) 

The estimated ratio of the staggered field for H||b and H|| c" was inconsistent between the specific heat and ESR. 
However, it should be recalled that each analysis has its own problem. In the specific heat analysis, an apparent 
gap structure, which is unrelated to the staggered field, was observed for H|| a". Whatever the origin of this gap 
structure, it is natural to expect similar contributions also for other field directions. Unfortunately, we do not know 
how to estimate these effects at present. On the other hand, there are also contribution to the ESR linewidth from 
other sources than the staggered field, and the subtraction causes some uncertainty, in addition to the estimate of the 
linewidth itself. 

The constraint from various experimental result on the DM vector is summarized in Fig. ( pi] ) . For the case of neutron 
scattering polarization we have included an estimated error bar from the polarization analysis of Ref. [ ^) . We have not 
attempted to estimate error bars in the other cases. We see that a candidate DM vector (D a >>, D c n) ~ (0.13, 0.02) J, 
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which satisfies the most reliable requirement of the cancellation for H|j a" , is roughly consistent with all the constraints 
except for the ratio of the gap between H||b and H|| c". 

This suggests that the observed ratio of the gap was wrong, or was affected by other factors than the staggered 
field. More experimental data are needed to draw a reliable conclusion. 

VIII. SUSCEPTIBILITY 

In the first subsection we consider the staggered susceptibility, resulting from the application of a staggered magnetic 
field, using the mapping onto the sine-Gordon model. In the next subsection we combine this with the standard uniform 
magnetization of the S=l/2 Heisenberg model to obtain the total physical susceptibility. 



A. Staggered Susceptibility 



We first discuss the susceptibility of the sine-Gordon model. Writing the sine-Gordon Lagrangian in the form: 

1 



£ = -d v 4>d v 4> + 2ficos(V2TT(/)), 

and adopting units where the velocity is set to 1, we define the sine-Gordon susceptibility as: 

d 2 F 



X. 



dfi 2 ' 



(8.1) 



(8.2) 



wliprp F is the free energy. The groundstate energy (i.e. the T = free energy) is expressed in terms of the gap 
as:l 



En = - 



4^/3' 



•3) 



This determines the T — susceptibility using the exact relationship between the coupling constant, /i and the mass 
A: 



A = n 2/3 A 



wherel 



A = A 2(2tt) 1 / 6 = 27T 1 / 6 
This gives the T = susceptibility: 

X(0) 

The high-temperature susceptibility is given by: 

rP r°° 



rr(3/4)- 


2/3 


rr(i/6)i 


|r(i/4)_ 




[r(2/3)J 



4.82764. 



A 3 
9VSA 



fP poo 

X^4 / dr dx < cosv / 27r0(r, x) cos V2%(f)(0, 0) > . 

Jo J -oc 



This Green's function is normalized to l/2r (at T = r = 0). At finite r and T we have: 

1 1 



Vr 2 + t 2 b 



sin ^(r + ix) sin ^(r — ix) 



Here (3 = 1/T. Thus, the susceptibility becomes: 

rP nT 



dx / dr- 

-oo JO 



1/2- 



r(i/4) 



sin ^ (r + ix) sin ^ (r — ix) 



1/2 T 



r(3/4) 



8.75376 
T ' 



(8.4) 



■5) 



(8.6) 



(8.7) 



•9) 
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(This result for the integral can be obtained, by analytic continuation, from the general results of Schulz.c3) An 
integral equation determining the sine-Gordon free energy at finite T was given in Ref . [ ^8| . We may determine the 
susceptibility by differentiating twice with respect to fi. Note however that this integral equation actually determines 
F(T) — Eq so we must add the zero-temperature part of the susceptibility, given above. The resulting susceptibility, 
%A is plotted versus T/A in Fig. (fl2|). As expected, it agrees quite well with high tempe ratu re res ult ~ 8.73576/T 



3> » 



down to T ~ A. It has a maximum at about T ~ 0.5A. The T = value is given by Eq. 

Up to a multiplicative factor and logarithmic corrections, the sine-Gordon susceptibility essentially gives the stag- 
gered susceptibility of the S=l/2 chain, i.e. its response to a staggered field: 



X>(T,h,H) 



d 2 F 



(8.10) 



In order to determine this factor and estimate the logarithmic corrections we first consider the T = staggered 
magnetization of the S=l/2 chain. 

We refer to the staggered magnetization as m s : 



< 5; >= 



(-l) j m a 



In the continuum limit, 



m s oc< trg<r x >oc< cos(27ri?0) > 



.11) 



.12) 



Since this operator has scaling dimension irR 2 , a standard RG scaling argument gives the scaling of the staggered 
magnetization with staggered field: 



J J \J 



7i\R a /(2— wBT) 



(8.13) 



for some function T>. For weak fields the exponent is approximately 1/3. In a similar way to our analysis of the gap in 
the previous section, by combining the exact results for the xxz model with an RG analysis of the marginal operator 
we may determine the scaling of magnetization with staggered field in the case of zero uniform field: 



In 



J 



1/3 



m s -> D (h/J)V 3 

and determine the behavior of D(H/J) in Eq. flS.13 ) for small field: 

V{H/J) -> DilniJ/H)} 1 / 3 . 



.14) 



.15) 



m s obeys a standard RG equation relating a change in the cut-off energy scale, E, to a change in the coupling constant, 
A: 



(8.16) 



Working to linear order in the marginal couplings, as before, we set (3 « and use Eq. (|J). Using Eq. ( gig ), and 
lowering the cut-off scale to the gap, A, this gives: 



F 







' 47T 


-1/4 




f 







for some constant, F. Using Eq. (5.11), and ( 5.33| )) this can be written: 



FA, 



1/2 



1/3 



47T 

71 



-1-1/3 



Ax (A) 



.17) 



(8.18) 



where Aq is defined in Eq. (5.5) and determined in Eq. ( 5.33| ). As in the previous sections we obtain the various 
formulae from the different asymptotic scaling of A_l(A) in the 3 cases: H = 0, H >> A and H = with exchange 
anisotropy e, defined in Eq. (5.21). Using Eq. (5.12) we obtain Eq. (3.14) with: 
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D = FAl'\ 



Using Eq. Q5.16D we obtain Eq. flS.13| ) with: 



H 



J 



D(-)=D Q {2(l-27rRi)}-^(- 



H 



-2(l-27ri? 2 )/3 



Using Eq. (5.14) we obtain Eq. (3.15) with: 



(8.19) 



(8.20) 



.21) 



Using Eq. ( 5.22 ), in the limit A/ J — > 0, in Eq. ( 8.17 ) we obtain an expression for m s in terms of A with exchange 
anisotropy: 



m s -> F(2e)- 1 / 4 (A/J) 



sl/2-e/4 



.22) 



We may determine the constant F, by comparing to the exact result of Ref. [ |12J)]. From Eq. ( 5.28 ): 

m s -> C(R) < exp[2niR4>] >, (8.23) 



with the exact formula for C(R) give in Eq. ( 5.2q ). This expectation value is given in terms of the soliton mass, A 
in [pj|, Eq. (15), with: 



(8.24) 



Inserting a power of the spin-wave velocity, v by dimensional analysis, and taking the limit of small e, we obtain: 



< e 



viR i_ (4/3)7rr(3/4) / r(2/3)r(5/6) V 

4^ J 



16sin(7r/3)r(l/4) 



- 3 / 2 /^\ 1 /2-c/4 



This gives a result consistent with Eq. ( 8.22 ) and determines the constant, F, to be: 

F = 2^A^£|Z| [r(2 /3)r(5/6)]- 3 / 2 = 2^/(3^), 



(8.25) 



(8.26) 



where the constant Aq is defined in Eq. ( |5.5| ) and i ts va lue is given in Eq. (5.33). Here we have used the exact identity 
r(l/6)r(5/6) = 2vr. Hence the amplitude of Eq. (|t|) is given by: 

D = FAl /2 = 2A 2 /(3V3tt). (8.27) 

We thus obtain the T — staggered susceptibility of the S=l/2 chain by differentiating m s with respect to h: 

2AI ln 1/2 (J/A) 



for H = 0. Comparing to Eqs. ( |3.5] ) and (^^), we see that: 

ln 1/2 (J/A) 



Xs (T = 0,A) = 



2(2tt) 3 / 2 



XSG (T = 0,A). 



(8.28) 



(8.29) 



The susceptibility for T >> A (and = 0^-,follows from Eqs. (|8.8| )-(8.9). Here we use the exact result for the 
T = correlation function of the S=l/2 chain:tj 



< S x (r)S x (0) >- 



(lnr 



,1/2 



(2tt) 3 / 2 7 



(8.30) 



This differs from the correlation function of the sine-Gordon model by a factor of (lnr) 1/,2 /2(27r) 3 / 2 . [Note the factor 
of 4 difference in the susceptibilities due to the factor of 2 in the interaction term of the sine-Gordon Lagrangian 
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of Eq. (3.1).] Upon going to finite T and Fourier transforming at zero frequency and wave- vector, we expect the 



logarithmic factor to become: In 1 / 2 (J/T). Thus: 



.277904 In 1 / 2 (J/T) Xn^jJ/T) 
Xs (T » A) = = 2(27r)3/2 X5g(T » A). (8.31) 



Comparing to Eq. flB.29 ) suggests the heuristic formula: 



v rr m_ ln 1/2 [J/max(r,A)] 

X s (i,Aj= 2(27r) 3 / 2 XSW,A). (8.32) 

XSG /[2(27r) 3 / 2 ] « X s is plotted in Fig. ©. 

Inclu ding a small uniform field, 7?, only makes unimportant changes in these formulas. At T — 0, the power of A 



in Eq. ( 8.28| ) changes by a small amount; the argument of the logarithm changes to J/H and the amplitude by a 
factor of V 1 - 73 . For T » A we must distinguish two regimes, depending on the relative magnitude of H and T. For 



A << T « H, the power of T changes. On the other hand, for A << H << T we expect to obtain Eq. (3.31 ) 



B. Physical Susceptibility 

Above we considered the staggered susceptibility resulting from (independent) staggered and uniform fields. To 
make any comparison with experiments we must take into account that the effective staggered field is proportional to 
the uniform field, 

h = cH, (8.33) 

where the constant of proportionality, c is strongly dependent on field direction. The physical susceptibility is 
conveniently obtained from its thermodynamic definition: 

d 2 F 

where F is the free energy. Writing F as a function of uniform and staggered fields, we m ust s et h = cH before 



taking the iJ-derivative. We may calculate the free energy in the rotated spin basis of Eq. (2.4), used throughout 
this paper. Noting that the first derivative of F with respect to H or h gives the uniform magnetization, m u and 
staggered magnetization m s respectively (in the rotated basis), we obtain: 



<9m„. 2 dm,, „ dm 



dH dh dH 

where we have used: 

<9m„ <9to„ 



*phy S = ^ +c l^+ 2c W' ( 8 - 35 ) 



dh dH 



dT(m s ( T )m u {T)), (8.36) 



We may ignore the dependence on h of the first term and use the standard result for the uniform susceptibility of the 
5 = 1/2 chain, Xu- For low fields and temperatures this gives: 



dm u ^ 1 _ 1 
~dH ~ Xu ~* 2^7 ~ 7rV 



X°u - 7T- = -n, (8-37) 



independent of field and temperature. The second term in Eq. ( 3.35 ) is larger than the third so we approximate: 

Xphys ~X° u + c 2 Xs(h,T), (8.38) 

where Xs is the staggered susceptibility discussed in the previous subsection. Thus, wiicn measuring the physical 
susceptibility of the present system, one actually probes also the staggered susceptibility 2 ^. 

While the first term, the standard result for the susceptibility of the S=l/2 chain, goes to a finite constant, at T 
an H — ► 0, the second term, resulting from the effective staggered field, is highly singular. At zero field it blows up 
at T — > as .278 In 1 / 2 (J/T)/T. Thus, although it is multiplied by the small constant, c 2 , it eventually dominates for 
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low enough T. For any finite field, the divergence of the second term is cut off, at essentially the gap energy, A (ft), at 
a value of approximately .229 In 1 / 2 (J/ A) /A, as shown in Fig. (|l2|). At low fields the behavior looks quite similar to 
a paramagnetic impurity contribution. It can be distinguished from that, however, by its very strong field-direction 
dependence. The effect is largest for the field in the c-direction when the parameter c 2 s» .01. 

The experimental susceptibilityLj of Cu benzoate shows very peculiar behavior at low fields and temperatures. As 
the temperature is lowered, in low fields, the susceptibility grows. This effect is highly direction dependent with the 
biggest effect occurring for fields in the c-direction. This effect is cut off by the application of a field. Qualitatively, the 
experimental results look similar to our predictions. However, at a more quantitative level there are many differences. 
The anomalous low field low T part of \ is much larger in the experiments than in the theory. For instance, at zero 



field, a temperature of IK w J/18, and c 2 = .01, the value that we find for H||c, the second term in Eq. (8.35) is 
still smaller than the fi rst by a factor of about 1/2. The anomalous contribution could be somewhat larger due to 
the third term in flB.35 ) and the logarithmic correction. However, the anomalous term in x c observed in Cu benzoate 



at this zero field and low T is about six times larger than the Xu- This is perhaps too large to be explained by our 
theory which is purely one-dimensional. 

Furthermore, the detailed experimental dependence on T and H is quite complex. The maximum susceptibility 
occurs at a field of about 30G. and T=1.5K. At lower fields or temperatures the susceptibility decreases somewhat. 
At low fields, two extra low temperature peaks are observed in the susceptibility as a function of T (in addition to the 
normal one for an S=l/2 chain at T ss .6J). While it is tempting to try to identify the higher-T peak with the peak in 
Xs shown in Fig. (|l2]), and the lower-T peak with the onset of Neel order, both susceptibility peaks are broad, unlike 
what would be expected from a phase transition. (Indeed, no evidence for magnetic order was found from neutron 
scattering,!] down to T=.8K.) A strong frequency dependence of the low T and H susceptibility was also observed. 

We expect that a proper theoretical description of the low field and temperature susceptibility of Cu benzoate will 
require the inclusion of inter-chain coupling effects. The interchain super-exchange paths look complicated and it 
is not even clear what is the sign of the interchain coupling. We remark that a perfectly one-dimensional S=l/2 
antifcrromagnet with SU(2) symmetry broken only by the DM interaction has a disordered groundstate. This follows 
from our spin redefinition in Sec. II which maps the system into an easy plane xxz model which is well-known to 
have a disordered groundstate with power-law correlations. We may think of the spins as fluctuating primarily in the 
plane perpendicular to D with a tendency for the cross-product of neighboring spins to be parallel to D. This system 
has a U(l) symmetry. Interchain coupling would normally be expected to produce long-range order with both an 
anti-ferromagnetic component and a possible perpendicular ferromagnetic component. For example, for Daz, 

< Sj >= -m s (-i) j x + m u y. (8.39) 

See Fig. ©. n 

The standard mean field treatment of the interchain interactions^ would suggest a critical temperature of order 
the interchain coupling (assuming that it is antifcrromagnetic) . On the other hand, when a magnetic field is applied 
(not parallel to D), the U(l) symmetry is broken and the phase transition should disappear. m s and m u now becomes 
non-zero even in the purely one-dimensional system and are no longer order parameters for spontaneous symmetry 
breaking. Thus the application of a magnetic field smoothes out the phase transition in this system. Surprisingly, 
neutron scattering experiments in zero field have failed to detect a Neel transition. 
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FIG. 1. Crystal structure of Cu Benzoate. Filled circles are Cu 2+ ions, connected atoms are benzoate group and grey circles 
represent H 2 molecules. Unit cell is shown as a frame, and arrows indicate crystal axes, part of the figure. 
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FIG. 2. Enlargement of crystal structure near a Cu (black spheres) chain with O atoms of H2O (dark spheres) and those of 
benzoate groups (light spheres) . Note that the oxygen octahedra have two different orientations on staggered Cu atoms. 
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FIG. 7. The RG flows in the presence of a magnetic field. The turn occurs at an energy scale of O(H). 
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.. Critical exponent 2irR 2 as a function of the applied field H. Exact solution is compared with the leading terms 
1 — l/[21og (Ho/H)] in the RG analysis. We show a good fit to the numerical solution of the Bethe ansatz integral 
equations with Ho = y'327r 3 /e. This choice of Ho is 4 times the value in Ref. [ |l6| , which gives worse fitting. 
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FIG. 9. Spin-wave velocity v as a function of the applied field H , determined from Bethe Ansatz integral equations. 
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FIG. 10. Magnetization m as a function of the applied field H, determined from Bethe Ansatz integral equations. 
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FIG. 11. Estimate of the DM vector from various experiments. Each constraint gives a set of allowed DM vectors as a curve 
in D a ii — D c ii plane. The constraint from neutron scattering polarization is drawn with an assumed error of ±5°. 
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FIG. 12. Susceptibility of the sine-Gordon model, as denned by Eq. |u|), (fj), divided by a factor of 2(2tt) 3/2 . This is 
essentially the staggered susceptibility of the spin chain, multiplied by A, up to a slowly varying logarithmic factor. The 
exact curve is obtained by a numerical solution of the integral equation, and the high-t emp e ratu re asymptotics is from the 
perturbation theory, .278A/T. The T = value is given by .229, in agreement with Eq. (O), (O). 
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